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H 要 : 全 球 气候 变 暖 背景 下 ,干旱 问题 更 加 频繁 持久 ,影响 范围 也 逐步 扩大 。 采 用 干旱 指数 进行 
干旱 评估 与 研究 时 , 受 多 种 因素 影响 ,干旱 评估 结果 往往 存在 一 定 的 不 确定 性 。 以 黑河 流域 为 研 
SIX. ,分析 概 率 分 布 模型 和 参数 估计 误差 这 两 种 不 确定 性 来 源 对 标准 化 降水 蒸 散 指数 (Standard- 
ized precipitation evapotranspiration index,SPFI) 和 于 时 特征 变量 (干旱 强度 .干旱 峰值 和 干旱 历时 ) 
的 影响 。 结 果 表 明 :这 两 种 来 源 均 会 对 SPEI 和 干旱 特征 变量 产生 影响 , 且 SPEI 越 极端 ,其 影响 越 
大 ,二 者 对 于 极端 和 严重 干旱 的 影响 程度 远大 于 对 轻 度 和 中 度 干 旱 的 影响 ;对 于 极端 和 严重 干旱 ， 
概率 分 布 模型 导致 的 干旱 强度 和 干旱 峰值 的 不 确定 性 更 大 ,参数 估计 误差 导致 的 干旱 历时 的 不 确 
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定性 更 大 。 研 究 结果 可 为 干旱 的 准确 评估 提供 支撑 ,在 防 旱 减 灾 工 作 中 为 制定 更 加 准确 有 效 的 搞 
旱 决 策 提 供 理论 支持 , 以 避免 可 能 造成 的 减灾 能 力 不 足 或 抗旱 资源 的 浪费 。 
X 键 词 : 概率 分 布 模型 ; 参数 估计 ; 不 确定 性 ; 标准 化 降水 蒸 散 指数 (SPEI) ; 黑河 


文章 编号 : 


气候 变 暖 是 人 类 面临 的 最 为 严重 的 全 球 环境 
问题 ,引起 了 世界 各 国政 府 和 公众 的 广泛 关注 "。 
全 球 变 暧 背景 下 ,干旱 问题 已 成 为 生态 学 、 气 象 学 
和 水 文学 领域 关注 的 重点 问题 之 一 。 干 旱 指数 是 
开展 干旱 研究 的 一 项 重要 工具 。 由 于 样本 、 模 型 以 
及 参数 的 影响 ,利用 干旱 指数 进行 干旱 的 监测 与 评 
估 时 ,得 到 的 结果 通常 具有 一 定 的 不 确定 性 。 

以 往 人 研究 中 ,样本 不 确定 性 对 于 干旱 指数 的 影 
响 研 究 相 对 较 多 。 例 如 , 洪 兴 骏 等 “采用 Monte- 
Carlo 随机 模拟 方法 验证 了 干旱 指数 (Streamflow 
drought index, SDI) 计 算 的 抽样 不 确定 性 ; 张 瑞 涵 等 
采用 参数 和 非 参 数 Bootstrap 方 法 评价 了 样本 不 确 
定性 对 标准 化 降水 指数 (Standardized precipitation 
index,SPI) 的 影响 。 在 模型 方面 ,干旱 指数 计算 时 
需要 用 到 概率 分 布 模型 来 拟 合 相关 数据 ,不 同 地 区 
气候 条 件 的 差异 导致 计算 时 所 选用 的 最 优 理论 概 
率 分 布 模型 各 不 相同 ;有 研究 表明 ,Pearson-II 分 布 
是 美国 地 区 以 及 我 国 孝 阳 湖 流域 计算 SPI 的 最 佳 模 
型 汪 ,欧洲 地 区 Weibull 分 布 比 Gamma 分 布 的 拟 合 
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效果 更 好 ;广义 正 态 分布 更 适合 巴西 地 区 多 时 间 
尺度 SPI 的 计算 ;欧洲 地 区 建议 使 用 广义 极 值 
(Generalized extreme value , GEV ) 分 布 计算 标准 化 降 
IK AS ACFE QL (Standardized precipitation evapotranspi- 
ration index, SPEL) 等 。 在 参数 方面 ,以 往 研究 更 多 
针对 不 同 参数 估计 方法 的 影响 ,Begueria 等 "对 比 了 
两 种 参数 估计 方法 一 一 最 大 似 然 估 计 法 和 无 偏 概 
率 加 权 算 法 在 计算 SPEI 时 的 性 能 。 而 模型 参数 全 
计 误 差 本 里 以 及 不 同 理论 概率 分 布 模型 对 干旱 评 
佑 结果 带 来 的 不 确定 性 在 现 有 研究 中 较 少 涉及 。 
鉴于 此 ,本 文 主要 探讨 不 同 概率 分 布 模型 以 及 参数 
估计 误差 这 两 种 不 确定 性 来 源 对 干旱 指数 以 及 干 
旱 评估 的 影响 程度 ; 由 于 SPEI 既 考虑 了 降水 、 又 考 
虑 了 蒸 散发 对 干旱 的 影响 ,因此 选择 该 指数 作为 待 
分 析 干 旱 指 数 。 

黑河 流域 是 我 国 西北 地 区 第 二 大 内 陆 河 流域 ， 
干旱 始终 是 该 流域 最 主要 的 自然 灾害 之 一 ”。 近 
2000 年 来 流域 累计 发 生 旱 涝 灾害 达 251 次 ,旱灾 占 
到 70% 以 上 "1!。 全 球 变 暧 背景 下 ,流域 中 下 游 的 干 
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旱 现 象 日 益 严 重 ,2001 年 发 生 了 1941 一 2001 年 近 
60 a 以 来 最 严重 的 秋冬 AEE, MAT 
早 灾害 不 仅 致使 当地 自然 生态 环境 异常 脆弱 ,也 进 
一 步 制 约 了 经 济 社会 的 健康 发 展 。 黑 河流 域 中 游 
走廊 平原 是 良好 的 发 展 农业 地 区 ,然而 由 于 受到 亚 
洲 东 岸 西北 气流 和 大 地 形 边缘 下 沉 气流 的 双重 控 
制 影响 ,该 地 区 是 全 球 同 纬度 有 人 类 活动 的 最 干旱 
的 地 区 之 一 ,农作物 需 水 高 峰 期 在 每 年 4 一 6 月 , 同 
时 期 正 是 该 流域 的 枯水期 ,导致 夏 灌 " 卡 脖子 " 旱 的 
矛盾 异常 突出 号 。 研 究 以 黑河 流域 为 例 , 评 估 了 9 
种 概率 分 布 模型 以 及 模型 参数 估计 误差 对 SPEI 以 
及 对 干旱 强度 干旱 峰值 干旱 历时 等 特征 变量 的 
影响 程度 ,从 而 为 干旱 的 准确 评估 提供 支撑 ,为 防 
旱 减 灾 工 作 中 更 加 准确 有 效 的 抗旱 决策 提供 理论 
支持 ,以 避免 可 能 造成 的 减灾 能 力 不 足 或 抗旱 资源 
的 浪费 。 


1 研究 区 概况 


黑河 流域 位 于 河西 走廊 中 部 ,地理 位 置 介 于 
97°37'~102°06’E, 37°44'~42°40’/N 之 间 ,流域 面积 近 
14.29x10' km ,流域 及 站 点 位 置 如 图 1 所 示 。 受 西 
风 带 环流 和 极地 冷气 团 的 双重 影响 ,研究 区 气候 存 
在 明显 的 南北 差异 与 东西 差异 ,跨越 了 3 个 地 理气 
候 区 ,由 南 向 北 分 别 为 : 半 王 旱 半 湿润 气候 区 SEE 
旱 气 候 区 和 极度 干旱 的 沙漠 区 ""。 上 游 南 部 山区 
多 年 平均 降水 量 400~500 mm, 多 年 平均 气温 1.5~ 
2.0 ,气候 冷 凉 ;中 游 多 年 平均 降水 量 100~150 
mm ,多 年 平均 气温 2.8~7.6 %C ;下游 多 年 平均 降水 量 
则 小 于 50 mm, 下 游 额济纳 平原 多 年 平均 气温 为 
8.0%。 总 的 来 看 ,全 流域 降水 稀少 但 分 布 集中 ,日 


À 


图 1 黑河 流域 及 气象 站 点 位 置 示意 图 
Fig. 1 Location diagram of the Heihe River Basin and the 


meteorological stations 


照 时 间 长 ,蒸发 作用 强烈 , 层 夜 温差 大 。 近 年 来 ,从 
流域 整体 来 看 ,其 平均 气温 呈现 出 显著 上 升 趋势 ， 
以 额济纳 旗 站 为 代表 ,其 增长 速率 达到 了 0.42 C- 
(10a)"” ,一 定 程度 上 加 剧 了 流域 的 干旱 问题 。 


2 数据 与 方法 


2.1 数据 来 源 

选用 流域 内 9 个 国家 级 气象 站 作为 研究 站 ,站 
点 基本 信息 见 表 1。 考 察 数据 的 完整 与 连续 性 后 ， 
选用 了 1960 一 2015 年 逐日 降水 量 和 气温 数据 集 用 
于 SPEI 的 计算 ,数据 来 源 于 国家 气象 科学 数据 中 心 


表 1 黑河 流域 气象 站 基本 信息 


Tab.1 Basic information of the meteorological stations in the Heihe River Basin 


气象 站 点 纬度 /N ZEE 1960—2015 年 多 年 平均 降水 量 /mm 1960—2015 EZE PEA a/c 
祁 连 站 38°11" 100?15' 409 INI 
7 牛 沟 站 38?25' 99?35' 418 -2.8 
TESI 38°48" 98°25" 298 -2.5 
山 丹 站 38°48" 101°05’ 201 6.6 
张掖 站 39°05’ 100?17' 128 7.5 
高 台 站 39°22" 99°50" 108 7.9 
酒泉 站 39*46' 98?29' 86 7.6 
鼎新 站 40°18" 99?3]' 55 8.5 
和 额济纳 旗 站 41°57’ 101°04" 35 9.0 
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(http://data.cma.en/) ,数据 质量 良好 。 

2.2 研究 方法 

2.2.1 45 E AUEK th ta AC (SPEL) do FF d EE 
€  Vicente-Serrano $E ^ f£ 2010 4E LA SPI 的 概念 和 
计算 方式 为 基础 ,提出 了 SPEI, 与 SPI 相 比 ,SPEI 考 
虑 了 温度 对 干旱 的 影响 ,常用 于 气温 敏感 地 区 的 干 


更 为 敏感 。 两 种 检验 的 计算 方法 详 见 Fadhilah 5 '* 

研究 。 

2.2.3 多 元 正 态 法 “采用 多 元 正 态 法 量化 由 参数 佑 

计 误 差 导 致 的 干旱 指数 估计 值 的 不 确定 性 "”””, 具 
(1) 使 用 LL3 概率 分 布 模型 分 别 拟 合 各 月 降水 


旱 特 征 分 析 , 并 建议 采用 三 参数 的 Log-Logistic 概率 
分 布 模型 来 描述 降水 量 与 潜在 燕 散 发 之 间 的 差 值 
变化 。 该 指数 的 计算 过 程 详 见 Vicente-Serrano 等 "9 
研究 。 根 据 SPEI 值 的 大 小 ,干旱 可 分 为 5 个 等 级 
(#2). 


#2 基于 SPEI 的 干旱 等 级 划分 
Tab.2 Drought classification based on standardized 
precipitation evapotranspiration index (SPET) 


干旱 等 级 SPEI 值 
无 早 —0.5«SPEI 
轻 度 干旱 —1.0<SPEI< -0.5 
中 度 干旱 —1.5«SPEI« -1.0 
严重 干旱 —2.0<SPEI< -1.5 
极端 干旱 SPEI< -2.0 

注 :SPEI 为 标准 化 降水 蒸 散 指数 。 


干旱 事件 的 发 生 指 的 是 干旱 指数 低 于 某 一 设 
定 的 国 值 ,例如 SPEI 值 介 于 -1.5~-2.0 之 间 时 ,表明 
一 次 严重 干旱 事件 发 生 了 。 根 据 游程 理论 ,干旱 
历时 即 干旱 指数 介 于 国 值 水 平 的 时 间 ;干旱 强度 定 
义 为 一 次 干旱 事件 中 的 干旱 指数 累加 值 与 干旱 历 
时 的 比值 ;干旱 峰值 定义 为 一 次 干旱 事件 中 的 干旱 
指数 最 小 值 。 

2.2.2 概率 分 布 模型 与 拟 合 优 度 检 验 ”为 评估 不 同 
概率 分 布 模型 对 SPEI 估 算 的 影响 ,选取 9 种 模型 进 
行 计算 ,具体 包括 三 参数 分 布 Log-Logistic (LL3)、 
Log-Normal (LN3) .GEV ,Pearson-III (PE3 ) #1 Weibull 
(WEI3) ,两 参数 分 布 Normal , Logistic , Chi-Square 和 
Cauchy。 各 概率 分 布 模型 的 概率 密度 函数 见 参 考 文 
献 L13]。 采 用 最 大 似 然 估计 法 对 不 同 分 布 模型 
进行 参数 估计 ,采用 Kolmogorov- Smirnov (K-S) fl 
Anderson-Darling(A-D) 两 种 检验 法 对 分 布 模型 进行 
拟 合 优 度 检验 。K-S 检 验 法 是 在 经 验 分 布 函数 的 基 
础 上 ,计算 给 定 累 积分 布 函数 和 经 验 分 布 函数 的 差 
值 最 大 值 来 进行 拟 合 优 度 检 验 , 这 种 方法 的 优点 是 
不 需要 事先 假设 拟 检验 数据 服从 的 具体 分 布 ;A-D 
检验 是 K-S 检 验 的 一 种 改进 , 它 对 分 布 两 端的 离 群 值 


ti (P) 5 E E Z& B (PET) B5 25 fe 38 FR y. 9] 
X = (s. Bat Xn) ,得 到 最 优 参数 集 (o. ^. y") 以 
及 对 应 的 方差 - 协 方差 矩阵 ,其 中 1 表示 月 份 ,m 表 示 
样本 量 大 小 。 

(2) 在 渐 近 正 态 假设 的 最 优 参数 集 和 方差 - 协 
方差 矩阵 的 基础 上 ,随机 生成 YV 组 参数 集 
os B. Sables Bos Ya ets yn) ,其 中 w 取 
1000。 

(3) 计算 每 组 参数 对 应 的 SPEI,V 组 参数 共 对 
应 N 组 SPEI。 

(4) 基于 N 组 SPEI, 定 义 第 (a/2) 和 第 (1-a/2) 个 
百 分 位 数 分 别 为 SPEI 的 置信 区 间 (CIs) 下 限 和 上 
限 ,“ 取 0.05 ,该 置信 区 间 即 是 考虑 参数 估计 误差 影 
响 的 SPEI 的 不 确定 性 区 间 。 


3 结果 与 分 析 


3.1 概率 分 布 模型 及 其 影响 
3.1.4 拟 合 优 度 检 验 结 果 首先 采用 Thornthwaite 
方法 计算 各 站 点 的 逐 月 潜在 蒸 散发 ,以 此 为 基 
础 ,使 用 不 同 概率 分 布 模型 拟 合 降水 量 与 潜在 蒸 散 
发 的 差 值 序列 。 图 2 以 1 月 为 例 , 给 出 了 不 同 概率 
分 布 模型 对 差 值 序列 的 拟 合 效果 。K-S 检 验 和 A-D 
检验 在 0.05 显著 性 水 平 下 的 临界 值 分 别 达 到 0.18 
£12.50 , 若 统计 值 低 于 这 一 临界 值 , 则 说 明 该 数据 序 
列 在 0.05 显著 性 水 平 下 通过 了 相应 检验 ,反之 则 表 
明 没 有 通过 。 

由 图 2 可 以 看 出 ,所 有 站 点 K-S 和 A-D 检 验 统计 
量 均 低 于 临界 值 , 说 明 这 些 分 布 都 可 以 用 来 拟 合 该 
序列 ,效果 良好 。 图 中 虚线 黑 框 表示 相应 站 点 统计 
量 的 最 小 值 ,可 以 看 出 ,9 个 中 有 8 个 最 优 模 型 是 三 
参数 概率 分 布 模型 ,其 中 LL3 拟 合 效果 最 好 ,对 应 的 
统计 量 最 小 值 最 多 , 占 到 61% ;其 次 是 LN3 和 WEI3 
分 布 。 两 参数 模型 中 ,Logistic 分 布 拟 合 效果 最 好 ， 
有 12 个 统计 值 最 小 , 占 比 达到 67%。Cauchy 分 布 的 
拟 合 效果 最 差 。 
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(a) K-S 检 验 


£X 0.08 
& 
$8 0.06 


147 (b) A-D 检 验 


=LL3 

a LN3 

a GEV 

a PE3 

= WEI3 

m Normal 

= Logistic 

= Chi-Square 
= Cauchy 


祁 连 站 ”野牛 沟 站 托 勒 站 “” 山 丹 站 ”张掖 站 ”高 台 站 ”酒泉 站 ” 易 新 站 额济纳 旗 站 
气象 站 点 


aLL3 

a LN3 

a GEV 

a PE3 

= WEI3 

m Normal 

1 Logistic 

= Chi-Square 
= Cauchy 


) 


祁 连 站 ”野牛 沟 站 托 勒 站 。 山 丹 站 ”张掖 站 ”高 台 站 酒泉 站 鼎新 站 额济纳 旗 站 
气象 站 点 


注 :各 站 点 K-S 和 A-D 检 验 的 统计 量 最 小 值 在 图 中 以 黑色 虚 框 表示 。LL3 、LN3 GEV .PE3 , WEI 分 别 为 
三 参数 的 Log-Logistic , Log-Normal , Generalized extreme value ,Pearson-III 和 Weibull 分 布 。 下 同 。 
图 2 不 同 概 率 分 布 模型 对 1 月 数据 序列 的 检验 统计 值 


Fig.2 Statistical values of tests of different probability distribution models for January data series 


此 外 ,其 余 月 份 数据 序列 的 拟 合 检 验 结果 也 表 

HH ,LL3 5 Logistic 分 布 分 别 是 三 参数 和 两 参数 模型 
中 拟 合 效果 最 优 的 分 布 模型 ,分别 有 48%(104) 和 
54%(116) 的 统计 量 最 小 。 
3.1.2 概率 分 布 模型 对 干旱 指数 的 影响 虽然 所 有 
分 布 均 适 用 于 拟 合 研究 区 降水 量 和 蒸 散发 的 差 值 
序列 ,但 不 同 分 布 计 算得 到 的 SPFI 是 否 接近 ,不 同 
分 布 的 选取 对 SPEI 带 来 的 不 确定 性 有 多 大 ,还 需要 
进一步 探究 。 由 此 采用 所 有 分 布 分 别 计算 各 站 点 
的 SPEI, 并 比较 各 分 布 与 最 优 分 布 (LL3) 计 算 结 果 
的 差异 。 

图 3 以 流域 上 游 时 和 牛 沟 站 ,流域 中 游 酒 泉 站 和 
流域 下 游 额 济 纳 旋 站 为 例 , 展 示 了 各 分 布 与 最 优 分 
布 计算 的 SPEI 的 残 差 值 散 点 图 。 可 以 看 出 ,SPEI 
等 于 0 为 轴 , 当 -1<SPEI<1 时 ,各 站 点 的 残 差 最 大 差 
异 值 介 于 -0.24~0.31 之 间 ; 随 着 SPEI 绝 对 值 的 增 
加 ,差异 也 随 之 增 大 ; 当 SPEI 增 加 至 +1.5( 严 重 干 


旱 ) 时 ,差异 显著 增加 ,最 大 残 差 值 高 达 -0.89~ 
-1.49。 与 其 他 分 布 相 比 ,Cauchy 导致 的 差异 最 大 
(黄色 曲线 ) ,Cauchy 分 布 计 算 的 指数 值 更 加 居中 ， 
极端 干旱 或 湿润 的 事件 较 少 ;例如 ,1992 年 3 月 野牛 
沟 站 和 1970 年 5 月 额济纳 旋 站 ,Cauchy 分布 计算 的 
数值 分 别 为 -1.37 和 1.45 ,属于 中 度 干 旱 .中 度 湿润 ， 
而 其 他 分 布 模型 计算 的 SPEI 的 平均 值 分 别 为 -2.09 
和 2.30, 属 于 极度 干旱 和 极端 湿润 。 男 外 , 除 Cauchy 
分 布 外 ,有 7 个 站 点 Logistics 分 布 给 出 的 差异 最 大 。 

图 4a 给 出 了 由 不 同 概率 分 布 模型 导致 的 所 有 
站 点 SPEI 的 差异 图 ,图 中 基本 呈现 出 随 着 SPEI 的 
绝对 值 的 增加 ( 越 湿润 或 者 越 干旱 ) ,差异 随 之 增 大 
的 趋势 ,尤其 对 于 中 度 干旱 (或 湿润 ) 以 上 等 级 。 当 
无 旱 和 轻 旱 时 ,SPEI 的 差异 范围 为 [0.03 ,0.45] ; 当 
中 度 干旱 (或 湿润 ) 时 ,差异 范围 为 [0.09,0.51]; 当 严 
重 干旱 (或 湿润 ) 时 ,差异 范围 增加 至 [0.26,0.96]; 极 
端 干 旱 ( 或 湿润 ) 时 差异 范围 为 [0.62 ,1.49]。 
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SPEI 指 数 
注 :SPEI 为 标准 化 降水 蒸 散 指数 。 
图 3 各 分 布 与 三 参数 Log-Logistic 分 布 计算 的 SPEI 的 残 差 值 
Fig. 3 Residuals of SPEI values derived from the alternative distributions and the 


optimal distribution (three-parameter Log-Logistic) 


FESD 。 野 牛 沟 站 。 祁 连 站 。 山 丹 站 。 张 掖 站 。 酒 泉 站 。 高 台 站 。 易 新 站 。 额 济 纳 旗 站 


1.6, 
L-6 [ (a) 不 同 概率 分 布 模型 对 SPEI 的 影响 : (6) 参数 估计 误差 对 SPEI 的 影响 
ê 
EN 1:2 
dal 
Ir 
ie 0.8 
tl 
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图 4 不 同 分 布 和 参数 估计 误差 对 SPEI 的 影响 


Fig.4 Effects of different probability distribution models and parameter estimation errors on the SPEI index 
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3.1.3 概率 分 布 模型 对 干旱 特征 变量 的 影响 基于 
游程 理论 提取 出 所 有 站 点 的 干旱 历时 .干旱 强度 和 
于 旱 峰 值 3 种 特征 变量 ,然后 按 表 2 的 干旱 等 级 划 
分 标准 进行 划分 。 图 5 给 出 了 不 同 干旱 等 级 下 所 有 
站 点 干旱 特征 的 统计 图 ,根据 统计 可 知 , 所 有 站 点 
轻 度 、 中 度 、 严 重 和 极端 干旱 强度 分 别 是 -0.74， 
-1.22,-1.69 和 -2.01。 由 不 同 分 布 导致 的 轻 度 .中 
BE 严重 和 极端 干旱 强度 的 93% 置 信 区 间 分 别 是 
[-0.81,-0.69] ,[-1.26,-1.10] .[-71.77,0]. [-2.36, 
0] , 随 着 干旱 等 级 的 增加 ,干旱 强度 的 不 确定 性 随 
之 增 大 ,尤其 对 于 严重 干旱 和 极端 干旱 ,其 95% 置 
信 区 间 宽 度 几 乎 达到 轻 度 和 中 度 干 旱 的 11~20 倍 
(图 5a)。4 个 等 级 干旱 峰值 的 95% 置 信 区 间 分 别 是 
[-1.00,-0.98], [-1.50, -1.34] .[ -2.00,0] ,[-2.48, 
0], 也 是 随 着 干旱 等 级 的 增加 ,其 不 确定 性 随 之 增 
大 ,严重 干旱 和 极端 干旱 峰值 的 不 确定 性 增加 尤为 
显著 (图 5b)。 对 于 干旱 历时 ,不 同等 级 干旱 程度 的 
不 确定 性 虽然 存在 差异 ,但 差异 不 大 , 轻 度 干旱 历 
时 不 确定 性 最 小 ,其 次 是 严重 干旱 和 极端 干旱 ,中 
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度 干 旱 历 时 不 确定 性 最 大 ,是 轻 度 干旱 不 确定 性 的 
2.7 倍 (图 $c)。 与 干旱 强度 和 王 旱 峰值 相 比 ,干旱 历 
时 的 不 确定 性 在 不 同 干旱 等 级 之 间 差 异 小 的 多 。 
3.2 参数 估计 误差 及 其 影响 
3.2.1 参数 估计 及 其 误差 ”利用 最 大 似 然 法 计算 得 
SI LL3 分 布 形状 参数 (shape) 尺度 参数 (scale) 和 位 
置 参 数 (location ) 参 数值 ,野牛 沟 站 TP 5g y RIED 
纳 旗 站 的 结果 如 表 3 所 示 , 表 中 以 1 月 为 例 ,分 别 给 
出 了 参数 最 优 值 和 参数 估计 值 的 95% 置 信 区 间 。 
3.2.2 参数 估计 误差 对 干旱 指数 的 影响 图 6 给 出 
了 3 个 站 点 SPEI 和 根据 多 元 正 态 法 得 到 的 SPEI 的 
95% 置 信 区 间 , 从 图 中 可 以 看 出 ,SPEI 越 大 或 越 小 ， 
对 应 的 置信 区 间 越 宽 。 

Fl 4b 给 出 了 由 参数 估计 不 确定 性 导致 的 SPFI 
的 95% 置 信 区 间 。 从 图 中 可 以 看 出 ,从 轻 度 干旱 
(或 湿润 ) 到 中 度 干 旱 ( 或 湿润 ) 再 到 严重 干旱 ( 湿 
润 ),95% 置 信 区 间 宽 度 明 显 呈 现 增加 趋势 。 当 无 旱 
和 轻 旱 时 ,95% 置 信 区 间 为 [0.53 ,0.79]; PREFE 
(或 湿润 ) 时 ,置信 区 间 为 [0.59,0.93] ;严重 干旱 (或 


mp BS% 
(a) 干旱 强度 (b) 干旱 峰值 
0.0 1 00r . 
-05 -0.5 
yc 
E -L0- 1} a -1.0 } 2s 
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干旱 等 级 干旱 等 级 
(c) 干旱 历时 
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注 : 黑 色 线 表示 由 不 同 分 布 导致 的 各 类 干旱 特征 变量 的 不 确定 性 ; 短 横 线 表示 干旱 特征 变量 的 最 值 ;原点 表示 均值 。 


图 5 不 同 分 布 和 参数 估计 误差 对 干旱 特征 变量 的 影响 


Fig.5 Effects of different probability distribution models and parameter estimation errors on drought characteristic variables 
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#3 三 参数 Log-Logistic 分 布 参数 最 优 值 及 9$% 置 信 区 间 


Tab.3 Optimal values and 95% confidence intervals for the parameters of three-parameter (Log-Logistic) distribution 


形状 参数 尺度 参数 位 置 参 数 
气象 站 点 95% 和 置信 区 间 95% 置 信 区 间 95% 置 信 区 间 
E 最 优 值 最 优 值 最 优 值 z 
下 限 EBR 下 限 上 限 下 限 ER 
野牛 沟 站 0.02 -0.18 0.21 34.79 27.11 42.37 68.96 51.86 85.52 
酒泉 站 -0.17 -0.39 0.03 25.50 19.65 31.70 —538.05 —550.86 -525.16 
额济纳 旗 站 0.04 -0.15 0.24 27.61 21.48 33.83 -704.63 -716.91 -691.00 
— SPEI  ---95%CIs 的 不 确定 性 ;对 于 干旱 历时 ,极端 .中 度 严重 干旱 


^T (a) 野牛 沟 站 


1960 1967 1973 1979 1985 1991 1997 2003 2009 2015 
年 从 


「 (b) 酒泉 站 


1960 1967 1973 1979 1985 1991 1997 2003 2009 2015 
年 份 


[额济纳 旗 站 


1960 1967 1973 1979 1985 1991 1997 2003 2009 2015 
年 份 


注 :CIs 为 置信 区 间 。 
图 6 SPEI 及 考虑 参数 估计 不 确定 性 的 95% 置 信 区 间 
Fig.6 Estimations of SPEI and the 95% confidence intervals 


with considering the uncertainty of parameter estimates 


湿润 ) 时 ,置信 区 间 增 加 至 [0.64,1.03]。 

32.3 参数 估计 误差 对 干旱 特征 变量 的 影响 图 5 
同时 给 出 了 由 参数 估计 不 确定 性 导致 的 不 同 干旱 
等 级 下 所 有 站 点 的 干旱 特征 。 轻 度 .中 度 .严重 、 极 
端 4 个 等 级 干旱 强度 的 95% 和 置信 区 间 分 别 为 
[-0.79, -0.72] .[- 1.30, - 1.14] .[ - 1.80, -0.66] 、 
[-2.28,0] ,干旱 峰值 的 置信 区 间 分 别 为 [-1.00， 
-0.98 ].[ -1.50,-1.40], | -2.00,-0.71].( -2.48,0], 
于 旱 历 时 的 置信 区 间 分 别 为 [7.43 ,13.15] ,[8.00, 
18.33] ,[1.28, 15.60] [0,15.73] ;干旱 等 级 越 高 , 干 
旱 强度 和 和 干旱 峰值 的 不 确定 性 越 大 ,极端 干旱 和 严 
重 干旱 的 不 确定 性 范围 明显 大 于 轻 度 和 中 度 干 旱 


历时 的 不 确定 性 范围 比较 接近 , 轻 度 干旱 历时 的 不 
确定 性 最 小 (图 5)。 


4 讨论 


根据 拟 合 优 度 检验 结果 ,选取 的 9 种 分 布 均 适 
用 于 研究 区 降水 量 与 潜在 蒸 散 发 差 值 序列 的 拟 合 ， 
三 参数 模型 中 LL3 分 布 适用 性 最 强 , 两 参数 模型 中 
Logistic 分 布 适用 性 最 强 ( 图 2); 如 果 仅 考虑 拟 合 效 
果 ,Logistic 分 布 由 于 其 仅 包 含 2 个 参数 ,计算 相对 
简单 ,因此 可 以 作为 拟 合 降 水 量 与 潜在 蒸 散发 差 值 
序列 的 候选 分 布 ; 但 通过 比较 SPEI 可 以 发 现 ,基于 
Logistic 分 布 与 基于 LL3 分 布 计算 得 到 的 指数 值 在 
极 值 处 存在 较 大 差异 ;虽然 两 参数 模型 参数 少 、 计 
算 简 单 ,但 该 模型 会 对 干旱 指数 的 计算 带 来 较 大 不 
确定 性 (图 3、 图 4)。 

通过 计算 不 同 概率 分 布 模型 对 应 的 SPEI 可 以 
发 现 , 概 率 分 布 模型 是 导致 SPEI 存 在 不 确定 性 的 来 
源 之 一 。Vergni 等 2 在 其 研究 中 也 发 现 , 除 数据 序 
列 长 度 外 ,概率 分 布 模型 是 导致 SPI 存 在 不 确定 性 
的 主要 来 源 ;Laimighofer 等 ”认为 ,概率 分 布 模型 对 
SPI 和 SPEI 不 确定 性 的 相对 贡献 程度 较 大 。 本 研究 
还 发 现 ,干旱 ( 或 湿润 ) 等 级 越 高 ,由 不 同 分布 模 型 
导致 的 SPEI 的 不 确定 性 越 大 ,尤其 对 于 严重 和 极端 
干旱 ,不 确定 性 增加 尤为 显著 ,甚至 会 影响 到 干旱 
等 级 的 评估 (图 3、 图 4)。 男 外 ,分 布 模型 还 对 干旱 
强度 .干旱 峰值 干旱 历时 等 特征 变量 的 评估 也 带 
来 一 定 的 不 确定 性 。 随 着 干旱 等 级 的 增加 ,由 其 导 
致 的 干旱 强度 .干旱 峰值 的 不 确定 性 均 随 之 增 大 ， 
严重 和 极端 干旱 的 不 确定 性 增加 尤为 明显 ,其 不 确 
定性 范围 达到 轻 度 和 中 度 干 时 的 10 倍 以 上 (图 5); 
这 说 明 由 概率 分 布 模型 导致 的 严重 和 极端 干旱 强 
度 和 干旱 峰值 的 不 确定 性 很 大 ,不 可 忽视 。 其 导致 
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的 干旱 历时 的 不 确定 性 并 未 随 干 旱 等 级 的 变化 而 
发 生 显著 变化 ,不 同等 级 之 间 干 旱 历 时 的 不 确定 性 
范围 比较 接近 (图 $)。 

模型 参数 估计 误差 也 会 导致 SPEI 存 在 一 定 的 
FEE”. Zhang 等 研究 表明 ,SPI 干 量 指数 
越 极端 ,由 模型 参数 估计 误差 导致 的 该 指数 不 确定 
性 越 大 ;本 文选 用 的 是 SPEI, 研 究 结论 与 其 类 似 , 即 
SPEI 越 极端 ,由 模型 参数 估计 误差 导致 的 SPEI 不 确 
定性 越 大 (图 4、 图 6)。 男 外 ,干旱 等 级 越 高 ,由 参数 
佑 计 误 差 导致 的 干旱 强度 .干旱 峰值 的 不 确定 性 也 
越 大 ,极端 干旱 和 严重 干旱 的 不 确定 性 范围 明显 大 
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参数 Logistic 分 布 都 能 较 好 地 拟 合 黑河 流域 的 
SPEI。 三 参数 概率 分 布 模型 的 拟 合 效 果 比 两 参数 
的 拟 合 效果 更 优 ,两 参数 模型 参数 少 、 计 算 简单 , 虽 
然 可 作为 拟 合 降水 量 与 潜在 蒸 散发 差 值 序列 的 备 
选 分 布 , 但 在 计算 SPEI 时 容易 引起 较 大 的 不 确定 
性 。 随 着 干旱 等 级 的 增加 ,由 不 同 概率 分 布 模型 导 
致 的 SPEI 的 不 确定 性 会 随 之 增 大 ,从 而 影响 干旱 评 
佑 结果 的 可 靠 性 。 

(2) 参数 估计 误差 也 会 导致 SPEI 及 干旱 评估 
结果 存在 一 定 的 不 确定 性 。SPEI 越 极端 ,模型 参数 
估计 误差 导致 的 SPEI 的 不 确定 性 越 大 ;干旱 等 级 越 


于 轻 度 和 中 度 干旱 的 不 确定 性 范围 ;不 同等 级 之 间 
干旱 历时 的 不 确定 性 范围 比较 接近 (图 5)。 

与 概率 分 布 模型 的 结果 相 比 ,尤其 对 于 轻 度 、 
中 度 .严重 干旱 ,由 参数 估计 误差 导致 的 SPEI 的 
95% 置 信 区 间 更 宽 , 即 不 确定 性 更 大 ; 由 参数 估计 误 
差 导 致 的 严重 干旱 强度 及 干旱 峰值 .中 度 干 旱 历时 
的 不 确定 性 偏 小 ,而 严重 和 极端 干旱 历时 的 不 确定 
性 则 有 所 偏 大 ,其 余 干 旱 特征 变量 的 不 确定 性 范围 
则 比较 接近 。 这 也 说 明 ,两 种 不 确定 性 来 源 对 干旱 
特征 评 佑 均 会 产生 影响 ,尤其 对 严重 和 极端 干旱 的 
强度 和 峰值 的 影响 较 大 。 因 此 ,在 进行 干旱 评估 
时 ,需要 重点 考虑 最 优 概率 分 布 模型 的 选取 以 及 最 
佳 参 数 估计 方法 ,从 而 降低 其 导致 的 干旱 评估 结 
的 不 确定 性 ,为 精准 干旱 评估 预测 预报 提供 支持 。 

需要 指出 的 是 , 除 文 中 讨论 的 两 种 不 确定 性 来 
源 外 ,干旱 评估 过 程 中 还 存在 其 他 不 确定 性 来 源 ， 
如 不 同 潜在 蔡 散 发 估算 方法 等 ;Aadhar 等 “研究 指 
出 ,潜在 落 散 发 估计 方法 对 干旱 评估 最 终结 果 的 不 
确定 性 贡献 很 大 ;限于 数据 资料 原因 ,文中 仅 采 用 
Thornthwaite 方 法 计算 各 站 点 湾 在 蒸 散 发 ,因此 也 会 
对 评估 结果 带 来 一 定 的 不 确定 性 ,这 方面 将 在 以 后 
的 研究 中 予以 综合 考虑 。 


5 结论 
本 文 以 黑河 流域 及 SPEI 为 例 ,探讨 了 9 种 概率 
分 布 模型 和 参数 估计 误差 这 两 种 不 确定 性 来 源 对 
干旱 指数 以 及 对 干旱 强度 .干旱 峰值 .干旱 历时 等 
特征 变量 不 确定 性 的 影响 ,主要 结论 如 下 : 
(1) 概率 分 布 模型 的 选择 对 干旱 评估 的 可 靠 性 
具有 重要 意义 。 总 体 而 言 ,三 参数 Log-Logistic 与 两 


高 ,干旱 强度 和 干旱 峰值 的 不 确定 性 越 大 ,极端 干 
旱 和 严重 干旱 的 不 确定 性 范围 远大 于 轻 度 和 中 度 
干旱 的 不 确定 性 。 

(3) 两 种 不 确定 性 来 源 对 干旱 评估 均 会 产生 影 
啊 ,尤其 对 严重 和 极端 干旱 的 强度 和 峰值 的 影响 较 
大 。 与 分 布 模型 相 比 ,参数 估计 误差 导致 的 轻 度 、 
中 度 .严重 干旱 的 SPEI 的 不 确定 性 更 大 ,严重 和 极 
端 干旱 历时 的 不 确定 性 更 大 。 这 说 明 ,在 进行 干旱 
精准 评估 时 ,以 上 两 种 不 确定 性 来 源 均 需 考 虑 ; 通 
过 选取 最 优 概率 分 布 模型 以 及 最 佳 参数 估计 方法 ， 
可 以 降低 干旱 评 佑 结果 的 不 确定 性 。 
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RGSS: 两 种 不 确定 性 来 源 对 干旱 指数 SPEI 及 干旱 评估 的 影响 


Effects of two uncertainty sources on drought index of SPEI 
and on drought assessment 


HAN Rucun', ZHANG Ying", LI Zhanling' 


(1. School of Water Resources and Environment, China University of Geosciences, Beijing 100083, China; 


2. Shouguang Water Conservancy Bureau, Shouguang 262700, Shandong, China) 


Abstract: Drought issues are becoming more frequent and persistent due to global warming, and the impact 
range is gradually expanding. Notably, drought monitoring, assessment, and research generally use drought 
indices. However, during their calculation processes, the sample, probability distribution function models, and 
parameters influence these indices leading to uncertainties in the drought assessment results. Therefore, this study 
investigated the influence of probability distribution function models and parameter estimation errors on the 
drought index- standardized precipitation evapotranspiration index (SPEI) and characteristics (drought intensity, 
peak, and duration) based on the SPEI in the Heihe River Basin in China as the study area. The results reveal that 
the two sources of the uncertainties, probability distribution function models, and parameter estimations, affected 
the drought index SPEI and assessment. Additionally, the results show that the more extreme the drought index 
SPEI is, the greater their effects, with both having a much greater impact on extreme and severe drought than on 
mild and moderate drought. The probability distribution function models resulted in more significant uncertainty 
in the drought intensity and peak of the extreme and severe droughts compared to the parameter estimation errors, 
which caused greater uncertainty in the drought duration in both drought conditions. Therefore, the study' s 
results can support the accurate assessment of drought and theoretically facilitate more accurate and effective 
decisions in drought prevention and mitigation efforts to avoid possible inadequacies in the mitigation capacity or 
waste of drought-resistant resources. 

Key words: probability distribution model; parameter estimation; uncertainty; standardized precipitation evapo- 
transpiration index (SPEI); Heihe River 


